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Large spin cold atom systems can exhibit novel magnetic properties which do not appear in usual 
spin-1/2 systems. We investigate the SU(i) resonating plaquette state in the three dimensional 
cubic optical lattice with spin-3/2 cold fermions. A novel gauge field formalism is constructed to 
describe the Rokhsar-Kivelson type of Hamiltonian and a duality transformation is used to study 
the phase diagram. Due to the proliferation of topological defects, the system is generally gapped 
for the whole phase diagram of the quantum model, which agrees with the recent numerical studies. 
A critical line is found for the classical plaquette system, which also corresponds to a quantum 
many-body wavefunction in a "plaquette liquid phase" . 



PACS numbers: 75.10 Jm, 75.40 Mg, 75.45.+j 



I. INTRODUCTION 

Quantum fluctuations and non-Neel ordering magnetic 
states in low dimensional spin-1/2 antiferromagnets are 
important topics in strongly correlated physics. The 
quantum dimer model (QDM) constructed by Rokhsar- 
Kivelson (RK) in which each dimer represents an SU (2) 
singlet provides a convenient way to investigate novel 
quantum magnetic states such as the exotic topologi- 
cal resonating valence bond (RVB) statesA. The QDM 
in the 2d square lattice generally exhibits crystalline or- 
dered phase except at the RK point where the ground 
state wavefunction is a superposition of all possible dimer 
coverings- 2 . In contrast, a spin liquid RVB phase has been 
shown in the triangular lattice in a finite range of inter- 
action parameters by Moessner et a! 3 . The 3d RVB type 
of spin liquid states have also been studied by using the 
QDMi£. 

Recently, there is a considerable interest on large spin 
magnetism with cold atoms in optical lattices 6 -' 7 -'^' 9 - : 1 ' 1 1 , 
whose physics is fundamentally different from its coun- 
terpart in solid state systems. In solid state systems, 
the large spin on each site is formed by electrons cou- 
pled by Hund's rule. The corresponding magnetism is 
dominated by the exchange of a single pair of spin-1/2 
electrons, and thus quantum fluctuations are suppressed 
by the large S effect. In contrast, it is a pair of large 
spin atoms that is exchanged in cold atom systems, thus 
quantum fluctuations can even be stronger than those in 
spin-1/2 systems. In particular, a hidden and generic 
Sp(A) symmetry has been proved in spin-3/2 systems 
without of fine-tuning by Wu et al££. This large sym- 
metry enhances quantum fluctuations and brings many 
novel magnetic physic o 9 ' 11 : 12 ' 13 : 14 . 

Below we will focus on a special case of spin-3/2 
fermions at the quarter-filling (one particle per site) in 
the 3D cubic lattice with an SU (4) symmetry which just 
means that all of the four spin components are equivalent 
to each other. The exchange model is the SU{A) antifcr- 
romagnetic Heisenberg model with each site in the fun- 
damental representation. Its key feature is that at least 



four-sites are required to form an SU (4) singlet two sites, 
i.e., two sites cannot form such a singlet. This SU(4) 
model was also constructed in spin-1/2 systems with or- 
bital degenerac y 15 : 16 . This model is different from the 
previous large- N version of the SU(N) Heisenberg model 
defined in the bipartite lattices where two neighbor- 
ing sites are with complex-conjugate representations and 
the Sp(2N) Heisenberg model defined in non-bipartite 
lattic o 17 ' 18 , both of which can have singlet dimers. The 
natural counterpart of the dimer here is the SU (4) singlet 
plaquette state as ^e a/ 3 7 sV4 (1)^(2)^7 (3)^1 (4), where 
a, (3, 7, and S take the value of S z as ±| and ±i. Re- 
cently, the crystalline ordered SU (4) plaquette state has 
been investigated in quasi-lD ladder and 2D square lat- 
tice system o 11 : 16 ' 19 . The resonating quantum plaquette 
model (QPM) in 3D has been constructed in Refi^S. where 
quantum Monte-Carlo simulation shows that the ground 
state is solid in the entire phase diagram. The SU(N) 
plaquette generalizations of the Affleck-Kennedy-Lieb- 
Tasaki (AKLT) states 21 have also been given in Refi22. 

In this article, we will formulate a novel gauge field 
representation to the resonating plaquette model based 
on SU (4) antiferromagnetic Heisenberg model in 3D cu- 
bic lattice. Unlike the QDM in 3d cubic lattice, this 
QPM is generally gapped for the whole phase diagram, 
due to the unavoidable proliferation of topological de- 
fects. We study the novel gauge field in the dual lan- 
guage, where a local description of topological defects is 
possible. The classical ensemble of the plaquette system 
is also discussed, and unlike its quantum version, our the- 
ory predicts the classical ensemble can have an algebraic 
liquid phase by tuning one parameter. Classification of 
topological sectors of the QPM is also discussed. 



II. QUANTUM PLAQUETTE MODEL 

The QPM model in the 3D cubic lattice can be rep- 
resented as follows. The effective Hilbert space is con- 
structed by all the plaquette configurations allowed by 
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the constraint: every site in the cubic lattice is connected 
to one and only one plaquette. Three flippable plaque- 
tte configurations exist in each unit cube corresponding 
to the pairs of faces of left and right, top and bottom, 
and front and back denoted as A, B and C in Fig. [TJ 
respectively. The RK-type Hamiltonian 1 reads: 

H = —t J2 {\A){B\ + \B)(C\ + \B)(C\ + h.c.} 
each cube 

+ V J2 {\A)(A\ + \B)(B\ + \C)(C\}, (1) 
each cube 

where t has been shown to be positive in Ref^S, and we 
leave the value of V/t arbitrary for generality. Eq. Q]can 
be represented as 

H = t ]T {|Qi)<Qi| + |Q 2 )<Q 2 |} + (^-2t) 
each cube 

x J2 {\A)(A\ + \B)(B + \C)(C\}, (2) 
each cube 

where |Qi) = \A) + u\B) + u 2 \C) , \Q 2 ) = \A)+lu 2 \B) + 
uj\C), and u — e l ^~ . As a result, at V — 2t (the RK 
point), the ground state wavefunction should be annihi- 
lated by the projectors |Qi)(Qi| and IQ2} {Q2 1 , i-e., the 
equal weight superposition between all the plaquette con- 
figurations which can be connected to each other through 
finite steps of local resonances, i.e., all the configurations 
within one topological sector. At V/t > 2 all the plaque- 
tte configurations without flippable cubes are eigenstates 
of the Hamiltonian, one of which is the staggered plaque- 
tte state. The phase diagram of this RK model has been 
studied numerically in Ref^i. In particular, the classi- 
cal Monte Carlo simulation performed shows that at this 
RK point a weak crystalline order of resonating cubes is 
formed which forms a cubic lattice with doubled lattice 
constant. At V/t < 2 the system starts to favor flippable 
cubes. For instance, at —V/t » 1 the ground states are 
twelve fold degenerate with columnar ordering. All the 
transitions between different phases are of the first order. 

The original RK Hamiltonian for quantum dimer 
model can be mapped to the compact U(l) gauge 
theory 2 * 2 ^, from which one can show that the 2+1 di- 
mensional QDM is gapped except for one special RK 
point, while 3+1 dimensional QDM has a deconfined al- 
gebraic liquid phase 24 . By contrast, the quantum pla- 
quette model in the cubic lattice can be mapped into a 
special type of lattice gauge field theory as follows. We 
denote all the square faces parallel to XZ plane of the 
cubic lattice by the sites to the left and bottom corner 
of the face: i + \x + and denote faces parallel to 
XY and YZ plane in a similar way. Then we define the 
boson number n with integer values on every face of the 
cubic lattice, n — 1 corresponds to a face with plaque- 
tte, and n = otherwise. A strong local potential term 
f7(fVl_i£-l_ij> — is turned on at every face to guaran- 
tee the low energy subspace of the boson Hilbert space is 
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FIG. 1: Three flippable configurations in one cube. The res- 
onance is represented in the t-terim in Eq. [Tj 

identical to the Hilbert space with all the plaquette con- 
figurations. Since every site is connected to one and only 
one plaquette, the summation of n over all twelve faces 
sharing one sites needs to be 1. Next, we define the rank- 
2 symmetric traceless tensor electric field on the lattice 
as 

E i>ilv = r/(i)(n i+ i A+ i e - i), Ep V = E VIM (n ^ u), (3) 

where = (— \y*+ l v+ l z equals 1 when i belongs to one 
of the two sublattices of the cubic lattice and equals — 1 
otherwise. It is straightforward to check that the one- 
site-one-plaquette local constraint on the Hilbert space 
can be written compactly as 

V x V y E xy + V y V z E yz + V Z V X E ZX = 5r]{i), (4) 

where V is lattice derivative with the usual definition 

V M / = /(i + A) -/(»)■ 

The canonical conjugate variable of Ei^ v is denoted as 
the vector potential of Ai^ u , 

Ai,^u = r)(i) fli+i/i+i*, Aftu = A vii , (n + 1 v). (5) 

^i+i/i+ip is the canonical conjugate variable of boson 
number n i+ i^ + i0, which is also the phase angle of boson 
creation operator. A^ and E^ satisfy 

[Ei,fj, u , Aj !p(T ] = i8i 3 (8pp8 ua + 8p a 8 vp ). (6) 

Because E^ v only takes values with an integer step, A^ 
is an compact field with period of 27r. Due to the com- 
mutator 

[Ei,fj, u , exp(iA, )I/(r )] = (SppSva + S^Sup) exp(iAj, va ), (7) 

operators exp(iA,- , va ) changes the eigenvalue of Ei^ v by 
1. As a result, the plaquette flipping process can be rep- 
resented as 

H t = - t{ cos{V z A xy ~ V x A yz ) (8) 
+ cos(V x A yz - V y A zx ) 
+ cos{V y A zx - V z A xy )}, (9) 
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which is invariant under the gauge transformation of 

A„ u -» A„ v + V M V„/, (10) 

which is already implied by the local constraint (|4|). / is 
an arbitrary scalar function. The low energy Hamiltonian 
of the system can be written as 



7?2 I rp2 \ 
cy 1 yz ' zx ) 



H = H t + U « 
each cube 

+ V {(^ x E yz ) 2 + (V y E zx ) 2 + (\7 Z E ZX ) 2 }, 

each cube 

(11) 

which is subject to the constraint in Eq. 2J Besides the 
gauge symmetry <jTUJ) , Hamiltonian (fTTj) together with 
constraint ^ share another symmetry as follows: 
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(14) 



/i, ^ and tr are three space coordinates. This symmetry 
forbids terms like E xy E yz to be generated under RG flow 
at low energy. 



III. DUALITY TRANSFORMATION 

A major question in which we are interested is whether 
the Hamiltonian Eq. [1] and Eq. [TT] have an intrinsic liq- 
uid phase, just like the 3d QDM in the cubic latticed. A 
liquid state here corresponds to a gapless Gaussian state 
in which we are allowed to expand the cosine functions 
in Eq. [TTJat their minima, i.e., a "spin wave" treatment. 
However, the Gaussian phase could also be a superfluid 
phase which breaks the conservation of boson numbers 
(or effectively the plaquette numbers) with (exp(i8)) ^ 0. 
In our current problem a superfluid phase is not possible 
because (exp(i6>)} ^ necessarily breaks the local gauge 
symmetry (ITU1) of Hamiltonian (fTTj) . In other words, a 
superfluid state is a coherent state of boson phase 6 im- 
plying a strong fluctuation of boson numbers, which obvi- 
ously violates the local one-site-one-plaquette constraint. 

In this type of lattice bosonic models, because bosonic 
phase variable A pl , is compact, the biggest obstacle of liq- 
uid phase is the proliferation of topological defect, which 
tunnels between two minima of the cosine function in Eq. 
IH1 Since the topological defects are nonlocal, the best 
way to study them is go to the dual picture, in which 
the topological defects become local vertex operators of 
the dual height variables, similar duality transformations 
have been used widely in studying all types of bosonic 
rotor models, such as in proving the intrinsic gap of 2D 
QDM 2,25 , showing the existence of "bose metal phase" 26 
as well as the deconfine phase of 3d QDM—, and very re- 
cently the stable liquid phase of three dimensional "gravi- 
ton" model22. 
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FIG. 2: The duality transformation defined in Eq. 1151 On 
dual site 1 and 2 there are three components of dual vector 
height /i M , and the dual transformation for the shaded face is 

Ey Z = (ky - h Z ) 2 - (ky ~ /i Z )l = V X {ky — ll Z ) . 



Besides the topological defects, another convenience 
one gains from the dual formalism is the solution of the 
constraint, i.e., we are no longer dealing with a Hilbert 
space with a strict one-site-one-plaquette constraint in 
Eq. UJ The dual variables are defined on the dual lattice 
sites i, which are the centers of the unit cubes. In order 
to solve the constraint completely, one needs to introduce 
three components of the height field (fi = 1, 2, 3) on 
every dual site i, which is the center of a unit cubic of 
the original lattice: 



E xy = V z (h x -h y )+E a xy , 
E yz = V x (h y -h z ) + E° z , 
E zx — V !/ (/i z — h x ) + E zx , 



(15) 



whose geometric illustration is shown in Fig. O h x , ViZ 
are fields only take discrete integer vaules. E xy ,E yz and 
E® x are background charges satisfying the constraint Eq. 
21 We can just take the configuration of the columnar 
phase to define the value of the background charges as 



-!!,.., J L-J— (when both i and j are even) 

~~ \ f_)i+i + fc + i , 

I v ' 2 (otherwise) 



E° xy (i,j,k) 



E°(i,j,k) = E° zx (i,j,k) = 



(-) 



i+j+k+l 



(16) 



The canonical momenta to the dual fields h p on 



each dual site are 



^ x V 7 y A zx ^7 z A xyj 
Tty ^ zA xy V x A yzj 
~Kz = ^ ' x A yz — W y A zx . 



(17) 



(18) 



One can check the commutation relation and see that 7r M 
and h p are a pair of conjugate variables. Then the dual 
Hamiltonian of (p~T|) reads 

H = - tcos ^ + Uj2C^ P [V»(h„-h p )~E„ p } 2 

fi—x,y,z fj,isp 

+ ^WV M (V M (A, - h p ) - El p )] 2 (19) 

fivp 

(,nvp is a fully symmetric rank-3 tensor which equals zero 
when any two of its three coordinates equal, and equals 
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one otherwise. On each dual lattice site i, the 7r M fields 
satisfy the relation that J2u= x y z k^j = 0. 

The symmetry transformations of Hamiltonian Eq. [19] 
can be extracted from the duality transformation Eq. [15] 
and Eq. QU 

h x ^h x + f(x, y, z) + gi(x), 
h y -> h y + f(x, y, z) + g 2 (y) 7 
h z ^h z + f(x,y,z)+g 3 (z), (20) 

where / is a function of three spatial coordinates and 
51,2,3 on ly depends on one spatial coordinate. This type 
of symmetry is a quasilocal symmetry, which also ex- 
ists in the Bose metal states 2 ^ and p-band cold atom 
systems 2 ^. 

The main purpose of this paper is to study whether 
Hamiltonian (|19p and (fTTj) have a liquid phase which pre- 
serves all the lattice symmetries, just like the deconfincd 
algebraic liquid phase of 3D QDM. In this kind of alge- 
braic liquid phase, one can expand the cosine functions 
in equation (|lip and relax the discrete values of the hp 
fields, the long distance physics can be described by a 
field theory which only involves coarsed grained mode of 
hp, let us denote the long scale mode as hp. In this Gaus- 
sian phase one can also define continuous tensor electric 
field Ep V as the coarse grained mode of Ep V , the rela- 
tion between Ep U and hp is Ep V = C,p yp d p {hp - h u ). 
A Gaussian field theory of hp should satisfy the con- 
tinuous version of symmetries listed in Eq. [^D] : hp — > 
hft + f{x,y,z) + .9^(77,), now hp as well as functions / 
and <? M can all take continuous values. A low energy field 
theory action is conjectured to be 

L = E 7 (W 2 + § E wv^, - h P )f + ■ ■ ■ , 

(21) 

where the h XtVyZ fields take continuous real values. No 
other quadratic terms of hp with second spatial derivative 
is allowed by symmetry in this action. Notice that in Eq. 
[5T]we have rescaled the space time coordinates to make 
the coefficients of the first and second term equal. The 
action (|2~Tj) describes a state with enlarged conservation 
laws of Tip. If there is a state described by the Gaussian 
action (|21|L 7r x , Tt y and -k z are conserved within each Y Z, 
ZX and XY plane respectively. So any operator with 
nonzero expectation values at this state has to satisfy 
the special 2d planar conservation law of Tfy. 

The Gaussian part of action (|2"Tj) has one unphysical 
pure gauge mode which corresponds to function / in Eq. 
1201 and two gapless physical modes, with low energy dis- 
persion : 

u\ ~ k 2 x + k 2 y + k 2 z + 

/i-4 T 1,4 T L4 _ U2 1,2 _ Ul 1,2 _ U2 1,2 
W K x -1- K y "I" K z K x K y KyK> z K X K Z 1 

«*l ~kl + k 2 + k 2 ~ 



l h.A _i_ £.4 _i_ ]U4 h.2 L.2 1.2 J.2 _ 1.2 1.2 (<yo\ 

y K x -t- K y ■+■ K z K x K y K y K z K X K Z . {/,/,) 

The second mode 0J2 vanishes at every coordinate axis 
of reciprocal space (k x ,k y ,k z ). The strong directional 
nature of uj 2 roots directly in the quasilocal gauge sym- 
metries in Eq. 1201 The same modes can be obtained from 
the continuum Gaussian limit action of Hamiltonian pip : 

L = E ^{(^A-) 2 - WM<^ - d„A p ,) 2 }. (23) 

In this action Ap V is the coarse grained mode of Ap V , 
and Ap V is no longer a compactified quantity. The fact 
that loi vanishes at every coordinate axis plays very im- 
portant role in our following analysis, since it will cre- 
ate infrared divergence along each axis in the momentum 
space, instead of only at the origin. Similar directional 
modes are also found in other systems with quasilocal 
symmetrie o 26 i 28 . 

The ellipses in Eq. [^contain the non-Gaussian vertex 
operators denoted as L v , which manifests the discrete na- 
ture of hp. Since hp only takes integer values, a periodic 
potential cos{2 / nhp) can be turned on in the dual lattice 
Hamiltonian (flT)]) . At low energy the Non-Gaussian term 
L v generated by cos(2nhp) has to satisfy all the symme- 
tries in Eq. [20] the simplest form it can take is cos[2nhp]. 
However, this vertex operator only has lattice scale cor- 
relation at the Gaussian fixed point, because it violates 
the gauge symmetry of action (|2ip . Thus the simplest 
vertex operator with possible long range correlation is 

L -» = E ~ acos [MK - K) + B^(i)]. (24) 

p=£v 

and B(i) is a function of dual sites, which is interpreted 
as the Berry's phase. The specific form of the Berry's 
phase of the vertex operators depends on the background 
charge of the original gauge field formalism, which de- 
termines the crystalline pattern of the gapped phase 25 . 
However, since the liquid phase is a phase in which the 
vertex operators are irrelevant, whether a liquid phase 
exists or not does not depend on the Berry's phase, thus 
in the current work we will not give a complete analysis 
of the Berry's phase of our problem. In the continuum 
limit the most relevant vertex operators are the ones with 
multi-defects processes without Berry's phase and consis- 
tent with symmetries cos[27riV(/L — h u )] 1 let us de- 
note this vertex operators as Vn„ v , and integer N can be 
determined from detailed analysis of the Berry's phase. 
The correlation function between two vertex operators 
with arbitrary N separated in space time is calculated as 
follows: 

(V N ^{0)V N ^(r)) 
~ exp{-(2^) 2 iV 2 ((/ lM (0) - M0)) (V(r) - h v (r)))} = 
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8 r ^r v exp{- 



2(27r) 2 iV 2 f d 4 k (2fc 2 + 3fc 2 +3fc 2 )e 



^ik-r 



K J (2tt)4 (Ag+ W ?)(*g+«2) J 

— > 8 r § Tu Const (r — > +oo](25) 



The correlation function (/i(r) h(r')) is evaluated at the 
Gaussian fixed point described by the continuum limit 
action (|21[) without L v . The delta function S riA S ril in 
Eq. [25] is due to the continuous quasilocal symmetry 
of action (f2Tj) . or in other words the conservation of 
7r M within each planes. For instance, correlation func- 
tion Ui 2 ^^{h x {0)-hy(0)} e -i^N{h x (r)-h y (r)}\ can omy ^ e 

nonzero when r x = r y = 0, otherwise tt x conservation 
within every YZ plane will be violated once r x ^ 0. 

Since the correlation function calculated in Eq. [231 
reaches a finite constant in the long distance limit, the 
vertex operators are very relevant at the Gaussian fixed 
point described by the action Eq. I21[ and the system 
is generally gapped with crystalline order in the whole 
phase diagram. Since this result is applicable to any N 
and independent of the Berry's phase, the same conclu- 
sion is applicable to all the QPM with a definite number 
of plaquette connected to each site. The specific crys- 
talline order can be determined from the detailed analysis 
of the Berry's phase. 



IV. CLASSICAL RK POINT 

At the RK point, the ground state wave-function is an 
equal weight superposition of all the configurations al- 
lowed by constraint Q . All the static physics of this state 
is mathematically equivalent to a classical ensemble, with 
partition function defined as summation of all the plaque- 
tte configurations with equal Boltzman weights. Since 
there is no energetic terms in the partition function, all 
that rules is the entropy. If we define the tensor electric 
field as Eq. [31 the classical ensemble can be written as 



(26) 
(27) 



The delta function enforces the constraint, and the term 
—U^ li _ LlJ (Ei^ v )' 2 in the exponential makes sure all the 
low energy E^ v configurations are one-to-one mapping of 
the plaquette configurations. Now solving the constraint 
by introducing dual height field the classical partition 
function can be rewritten as 

Z = 5>xp{-E/^£ WV^i,, - h, p ) - <] 2 X28) 

hj^ i fiup 

Again we are mainly interested in whether this classi- 
cal ensemble is an algebraic liquid state, or by tuning pa- 
rameters one can reach an algebraic liquid phase. We can 



conjecture a low energy classical field theory generated by 
entropy, allowed by symmetry ()20p . The same strategy 
has been used to study classical six-vertex model, classi- 
cal three-color model and four color model 2 ^. Here the 
simplest low energy effective classical field theory reads 



(29) 



The number K cannot be determined from our field the- 
ory. This is the simplest free energy allowed by sym- 
metry. The physical meaning of this free energy is that, 
the total number of plaquette configurations (entropy) in 
a three dimension volume is larger if the average tensor 
electric field E^ v is small, i.e. the entropy favors zero 
average tensor electric field. 

The ellipses in equation ([2"9"|) includes the vertex opera- 
tors in equation (|24[) . The relevance of the vertex opera- 
tors can be checked by calculating the scaling dimensions 
of the vertex operators at the Gaussian fixed point action 
([2"9"| . Let us denote vertex operator cos[27rA(ft, Al — h u )] 
as Vn.^v Due to the symmetry ([20)) , Vn, xv can only cor- 
relates with itself along the same z axis, and Vm,zx and 
Vn, V z can never have nonzero correlation between each 
other when they are separated spatially along z axis. 

The leading order correlation functions are 



.xy, (0,0,0) 



exp{- 



(2tt) 2 A 2 2 f d 3 k (kl + k 2 y y k * z 



K 



Ull,2 i 1,2 1,2 I 1,21,2 
1 



} 



z8 ^/*- (30) 

In the above calculations we have chosen the simplest 
rcgularization: replacing spatial derivative on the lattice 
by momentum ik x . It has been shown that the scaling 
dimensions of operators in these type of models with ex- 
treme anisotropy can depend on the regularization on the 
lattice^. Here the scaling dimension of operator Vn.^v is 
regularization independent. These vertex operators are 
irrelevant if K < K c = 4irN 2 , in this parameter regime 
the contribution of Vff }X y to various correlation functions 
can be calculated perturbatively. 

Some other vertex operators can be generated un- 
der renormalization group flow, but these vertex oper- 
ators all have algebraic correlations, with a regulariza- 
tion dependent scaling dimension proportional to 1/K. 

For instance, vertex operator cos[27rA f VJ^ a .(/i a ; — h y )] has 
nonzero algebraic correlation function in the YZ plane 
at long distance. Here lattice derivative V mx is de- 
fined as V mx f(r) = f(f+ mx) — f(r). If we regular- 
ize the theory by replacing lattice derivative V mx with 
i2 sm(mk x /2) in the momentum space, the scaling dimen- 
sion of VN,mx,xy = cos^TriWJJ^. (h x — h y )] with n = 1 and 
arbitrary integer m is 4:TrN 2 /K, and the scaling dimen- 
sion is isotropic in the whole YZ plane: 

c{6) 



,mx,xy, (0,0,0) 



fy2 _|_ z 2^4ttN 2 /K 



.(31) 



6 



Here c(6) is a positive function of 9 = arctan(z/y). No- 
tice that the rotation symmetry in the Y Z plane is not 
restored even at long length scale. The scaling dimension 
of cos[27riVV^ ix (/ia ; — h y )] increases rapidly with number 
n. Thus all the vertex operators are irrelevant when K 
is small enough, and there is a critical K c separating a 
crystalline order and the algebraic liquid phase. At the 
liquid line the crystalline order parameter should have al- 
gebraic correlation functions. Coefficient K can be tuned 
from adding energetic terms in the system. Recall that 
now the configurations with zero average E^ v are favored 
by entropy, if we want to reduce K, we can add energetic 
terms which disfavor zero average E^ v . For instance, if 
we give the flippable cubics a smaller weight than the 
unflippable cubics, coefficient K should be reduced. 

The above results can be roughly understood from sim- 
ple physical argument. Notice that all the flippable cubes 
have zero average electric field, so the entropy effectively 
favors flippable cubes. If K > K c , the entropy strongly 
favors flippable cubes, the system will develop crystalline 
order which maximize the number of flippable cubes. 
This kind of effect is usually called "order by disorder" . It 
is also natural that the crystalline order tends to be weak- 
ened or even melt if we reduce K. Since the melting tran- 
sition of the crystalline order is driven by the proliferation 
of defect operators, the universality class of this transi- 
tion is very similar to the Kosterlitz-Thouless transition 
of the 2D XY model, the correlation length between oper- 
ators Vn,hu diverges as £ ~ exp[^/ b/ {K — K c )]. Unusual 
KT like transition in 3d or higher dimensions have also 
been discussed in other systems with similar quasilocal 
symmetries^, where the dimensionality of the system is 
effectively reduced to 2d. 

Recent Monte Carlo simulation^ shows that the whole 
phase diagram of RK Hamiltonian |T]) is gapped with 
crystalline order, including the RK point. Our results 
based on duality is consistent with this numerical results, 
and the equal weight classical partition function should 
have K > K c . Our theory also predicts that if we turn 
on energetic terms which favors unflippable cubes, there 
is a critical line described by Gaussian field theory (f2U)) . 
This prediction can be checked by classical Monte Carlo 
simulations. Another prediction which in principle can 
be made in our formalism is the most favored crystalline 
order when K is slightly larger than K c . This requires a 
detailed analysis of the Berry's phase of the vertex oper- 
ators in the dual theory, which we leave to future studies. 



V. TOPOLOGICAL SECTOR 

Now let us discuss the topological sector, within which 
every configuration can be connected to each other 
through finite local movings depicted in Fig. [TJ Topolog- 
ical sectors are especially useful when one is dealing with 
a quantum liquid state, where Landau's classification of 
phases are no longer applicable. In the original quantum 



z 




a b 

FIG. 3: (a), the conserved quantity m xy is defined as the 
summation of all E xy on all the shaded squares along one 
2 axis. (6), the view of the 3d lattice from the top. If the 
quantity m xy is fixed on all the shaded squares shown in this 
figure, m xy is determined on the whole lattice. 

dimer model on square lattice, the topological sector on 
a torus is specified by two integers 1 , which can be in- 
terpreted as winding numbers of electric fields. Here we 
choose a lattice with even number of sites in each axis 
and impose the periodic boundary condition. To spec- 
ify a topological sector one needs to know the conserved 
quantities under local movings. It is straightforward to 
check that quantity TO^j xy = E ixy for any 2d co- 
ordinate (i x ,i y ) is a conserved quantity. Notation Y^ 
means summation over all the sites with the same x and 
y coordinates (Fig. [3]). However, these quantities are not 
independent. For instance, using constraint Q we have 
the following identity: 

mofi,xy — mi,0,xy + mi,l,xy — fn^X^xy = VxVyExy 

= 5r >W - & y V z E yz + VjVx-E,*) = 0.(32) 

i z 

Thus as long as one fix the quantity m xy for one column 
and one row in the XY plane, their values for the whole 
lattice are determined. Conserved quantities associated 
with E zx and E yz can be treated in the same way. Thus 
we conclude that one needs infinite number of integers to 
specify a topological sector on a three dimension torus, 
and the number scales with the linear size of the lattice. 



VI. SUMMARY AND COMPARISON WITH 
OTHER MODELS 

This work studies a three dimensional quantum res- 
onating plaquette model, motivated from a special SU(4) 
invariant point in spin-3/2 cold atom system. The effec- 
tive low energy physics of the problem can be mapped to 
a special type of lattice gauge field. Our current QPM 
together with previously studied 3d QDM21 and soft- 
graviton model22, all have local constraint and low energy 
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gauge field description without gapless matter fields. Un- 
like the QDM and the soft-graviton model, the QPM al- 
most always suffers from the proliferation of topological 
defects, and a generic stable algebraic liquid state as an 
analogue of the photon phase of 3d QDM does not exist. 

The reason of the existence of a stable liquid phase 
of 3d QDM as well as the 3d soft-graviton model have 
been discussed in referenced. Both models with stable 
liquid phases are self-dual gauge theories, with strong 
enough gauge symmetries in both the original description 
of the problem or the dual theories, i.e. one cannot write 
down a gauge invariant vertex operator which gaps out 
the liquid phase. In our current QPM, the symmetry of 
the dual theory does not rule out all the vertex operators, 
and gauge invariant vertex operators are very relevant. 



Thus in this type of bosonic quantum rotor models, large 
enough gauge symmetries are necessary for both sides of 
the duality to guarantee the existence of a stable liquid 
phase if gapless matter field is absent. 
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